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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


RESEARCH MEMORANDUM 


EFFECT OF WING HEIGHT AND DIHEDRAL ON THE LATERAL STABILITY 
CHARACTERISTICS AT LOW LIFT OF A 45° SWEPT-WING AIRPLANE 
CONFIGURATION AS OBTAINED FROM TIME-VECTOR ANALYSES 
OF ROCKET-PROPELLED-MODEL FLIGHTS AT 
MACH NUMBERS FROM 0.7 TO 1.5 
By Clarence L. Gillis and Rowe Chapman, Jr. 

SUMMARY 


Lateral-stability flight tests were made over the Mach number range 
from 0.7 to 1.3 of models of three airplane configurations having 45° 
sweptback wings. One model had a high wing; one, a low wing; and one, 
a high wing with cathedral. The models were otherwise identical. The 
lateral oscillations of the models resulting from intermittent yawing 
disturbances were interpreted in terms of full-scale airplane flying 
qualities and were further analyzed by the time-vector method to obtain 
values of the lateral stability derivatives. 

The effects of changes in wing height on the static sideslip deriv- 
atives were fairly constant in the speed range investigated and agreed 
well with estimated values based on subsonic wind-tunnel tests. Effects 
of geometric dihedral on the rolling moment due to sideslip agreed well 
with theoretical and other experimental results and with a theoretical 
relation involving the damping in roll. The damping in roll, when com- 
pared with theoretical and other experimental results, showed good agree- 
ment at supersonic speeds but was somewhat higher at a Mach number of 1.0 
and at subsonic speeds. The damping in yaw showed no large changes in 
the transonic region. 

The damping of the lateral oscillations for typical full-scale 
airplane flight conditions was marginal or unsatisfactory. In addition, 
the high-wing model had an undesirably high ratio of rolling motion to 
yawing motion. 
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INTRODUCTION 


Attainment of flight speeds in the supersonic realm has focused 
attention on the stability problems associated with flight at high speeds 
and high altitudes. The types of lateral motions encountered and the 
importance of certain of the stability derivatives and inertia relations 
in defining the motions of the airplane have been examined. (See refs. 1 
to 5, for example.) Techniques for determining the magnitude of the 
pertinent derivatives have been established for some time in the subsonic 
speed range. For supersonic speeds most of the tests have been made to 
determine the static lateral force and moment and the rolling moment due 
to sideslip and roll velocity. Recent work in the Langley 9 -inch super- 
sonic tunnel produced measurements of the rolling moment due to the rate 
of change of yaw and sideslip (ref. 6) . Other techniques and measurements 
for the damping and cross derivatives are given in references 7 and 8. 

In order to furnish additional information on the lateral stability 
derivatives at transonic speeds, a rocket-propelled-model program to deter- 
mine the effects of wing height and dihedral on the lateral stability 
characteristics of a 45 ° swept-wing airplane configuration has been con- 
ducted. The present paper reports the results obtained from three free- 
flight rocket-propelled models of this test program. - The data were 
obtained on a high-wing model, a low-wing model, and a high-wing model 
with 10° cathedral. 

An additional purpose of the present program was to provide further 
experience for evaluating the test technique and the time-vector method 
of data analysis. This method has been used previously for specific cases 
of both full-scale and model lateral-stability flight tests (refs. 9 to 12) . 
The geometric variations utilized in the program to be discussed herein 
provided appreciable variations in frequency and very large variations in 
roll-to-yaw amplitude ratios. Some evaluation of the effects of these 
factors on the usefulness and accuracy of the test technique was possible. 

The models were flight tested at the Langley Pilotless Aircraft 
Research Station at Wallops Island, Va. 
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SYMBOLS 


chord-force coefficient, 
lateral-force coefficient, 
normal-force coefficient. 
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rolling-moment coefficient. 


Rolling moment 
qSb 


n 

C L 

I v 


X xz 

% 

L 

M 


yawing-moment coefficient. Y awing moment 

<jbD 

lift coefficient, C-^ cos a - Cc sin a, 

p 

moment of inertia about X-axis, slug-ft 

O 

moment of inertia about Z-axls, slug-ft 

product of inertia, i(l z - lx) i an slug-ft^ 

span-load distribution factor, ft^ (see appendix B) 

concentrated load, lb (see appendix B) 

Mach number 


Mjj, My, M^ rolling, pitching, and yawing moments about X, Y, and Z axes, 
respectively 


P 

R 

S 

T l/2 

V 

W 

X, Y, Z 
a 


period of lateral oscillation, sec 

Reynolds number based on wing mean aerodynamic chord 
wing area, sq ft 

time required for lateral oscillation to damp to half 
amplitude, sec 

velocity, ft /sec 

weight of model, lb 

coordinate axes 


damping factor, 


-0.693 1 


T l/2 ' 


sec 


m 


longitudinal accelerometer reading, positive in the positive 
x-direction, g units 

normal accelerometer reading, positive in the negative 
z-direction, g units 
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b 

c 

c 

c l/2 


g 

*x 

^2 

l 

m 

P 

q 

r 

t 

v e 

y 

y* 

a 

3 


lateral accelerometer reading, positive in the positive 
y-direction, g units 

wing span, ft 

local chord, ft 

mean aerodynamic chord, ft 

number of cycles required for oscillation to damp to half 
amplitude 

gravitational acceleration, 32.2 ft/sec^ 
radius of gyration about X-axis, ft 

radius of gyration about Z-axis, ft 

rolling moment on a section of the wing, ft-lb (see appendix B) 
longitudinal distance between two transverse accelerometers, ft 
mass, slugs 

rolling velocity, radians /sec 

dynamic pressure, lb/sq ft, or pitching velocity, radians/sec 
yawing velocity, radians /sec 
time, sec 

equi valent lateral velocity, ft/sec 
spanwise station, ft 

spanwise station at which concentrated load is applied, ft 
angle of attack, radians or deg 
angle of sideslip, radians or deg 

dihedral angle, deg, or span- load parameter (see appendix B) 


r 
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V’b 


* 

^c Y 


co 

a*i 


principal-axis inclination, positive as shown in figure 5 


pitch angle or angle of twist (see fig. 3 a nd appendix B) , 
radians 

relative mass density, -Sj- 

pSb 

atmospheric density, slugs /cu ft 

roll angle, radians 

yaw angle, radians 

phase angle "by which Cy leads 0 

phase angle by which p leads 0 

oscillation frequency, radians/sec 
undamped natural frequency, radians /sec 


Subscripts : 

T trim value 

A full-scale airplane 

r rigid 

e elastic 

WF wing-fuselage 

V vertical tail 

W wing 

H horizontal tail 

A single dot over a symbol indicates the derivative of the quantity 
with respect to time; a double dot represents the second derivative with 
respect to time. Amplitude ratios of the oscillatory components of the 

motion are designated as |~L |^| , and so forth. The static stability 
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dCy 

derivatives are indicated in the following manner: Cy = — nnri so 

P 

forth, whereas the rotary and acceleration derivatives are indicated 
dCn Sen 

as : On,, = — Cns = and so forth. 

°2V ^ 


MODELS AND INSTRUMENTATION 
Models 


Drawings and dimensions of the models are presented in figure 1. 
Additional geometric characteristics of the models are given in table I 
ayiri the mass and inertia characteristics are given in table II. The 
three models had geometrically identical surfaces and fuselages and dif- 
fered only in the vertical position and dihedral of the wing. The wing 
weight was counterbalanced so that the model center of gravity was on the 
geometric center line for all models. 

Wings and tail surfaces had the quarter-chord line swept back 45° 
and had NACA 65A004 airfoil sections. The wings were composite wood 
and metal construction. Verical- and horizontal- tail surfaces were 
machined from solid duralumin. The horizontal-tail surface was set at 
an incidence angle of -l/ 2 ° with the fuselage center line for all models. 

The fuselage was a parabolic body of revolution described by the 
coordinates given in table III. Housed in the spun-metal nose of the 
model was an electrically actuated yaw vane disturber. Photographs of 
the models are shown in figure 2 . 

In order to minimize the amount of pitch-yaw coupling encountered 
in the tests, the center of gravity of each model was placed very far 
forward (see table II) in an attempt to keep the pitch frequency always 
higher than twice the yaw frequency and thus avoid a resonant condition 
(ref. 13 ) . 

Throughout the remainder of this report, model 1 will designate 
the high-wing model, model 2 , the low-wing model, and model 3 > the 
high-wing model with cathedral. 


Instrumentation 

The models each contained an eight-channel MCA telemeter. The 
quantities recorded continuously were total pressure, angle of sideslip, 
rolling velocity, normal, longitudinal, and transverse accelerations 
near the center of gravity, transverse acceleration at the nose. 
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yaw- vane -disturber deflection (model 1 only) , and angle of attack (models 2 
and 3 only) . 

Ground-based instrumentation included a CW Doppler radar unit for 
measuring model velocity and acceleration, a modified SCR -584 radar unit 
for obtaining the model position in space as a function of time, and a 
rollsonde receiver operating on the telemeter-antenna radiation pattern 
as an additional measure of the model rolling velocity. Atmospheric 
conditions and wind velocities over the firing range were obtained from 
rawinsondes released immediately after the model flights. Manually 
operated l6-millimeter and ^-mlllhaeter tracking cameras were used to 
photograph the model flight. 


TESTS AND ANALYSIS 
Preflight Tests 


Because of the thin airfoil sections used, the model lifting surfaces 
were quite flexible. In order to obtain stiffness data, one wing panel 
and the vertical tail of one model were subjected to concentrated static 
loads at four points along both the 0.25 and 0.50 chord lines and the 
deflections along the leading and trailing edges measured. The results 
of the static loading tests are presented as influence coefficients In 
figure 3. The moments of inertia were determined by swinging the model 
as a pendulum and measuring the frequency of the oscillations. The 
principal longitudinal axis was determined by swinging the model in roll 
about a number of longitudinal axes inclined at various angles in the 
plane of symmetry and noting the angle for which the moment of inertia in 
roll was a minimum. 


Flight Test 

Each model was accelerated to a maximum Mach number of approx:* 
imately 1.5 by a solid-fuel ABL Deacon rocket motor which separated from 
the model after burnout as a result of the relative drag-to-weight ratios. 
Figure 2 (c) shows the high-wing model and booster on the launcher. The 
models contained no sustainer rockets. The yaw disturbers remained 
retracted during boosted flight and their operation was started at sep- 
aration. The yaw vanes were programed to extend fully in a time of 
0.33 second and retract in a time of approximately 0.05 second, the cycle 
being repeated every 1.1 seconds. Since the yaw vanes were set at a fixed 
angle of 10° with respect to the model center line, the model assumed a 
negative sideslip angle with the vanes extended. Upon the rapid retrac- 
tion of the vanes, the model oscillated about the undisturbed trim 
conditions. The test conditions prevailing for the flights are pre- 
sented in figure 1 as a function of the Mach number. 
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Corrections 

Position corrections were required for some of the accelerometer 
readings to convert the measured values at the instrument location to 
values at the center of gravity. Corrections for model pitching and 
yawing velocities were also made to the readings of the air-flow indi- 
cators to obtain angles of attack and angles of sideslip. In the process 
of data analysis, all the measured lateral-oscillation data were refer- 
enced to 3 . The frequency-response characteristics of all the instru- 
ments were sufficiently close to those of the sideslip vane so that the 
only correction required was to the phase angle of the rolling velocity. 
This correction varied between approximately 1° and 5° for the tests 
herein. 


Analysis 

The flight-test data presented herein were analyzed by the time- 
vector method of solution of the lateral stability equations of motion. 

The equations of motion used in the analysis and the assumptions involved 
are given in appendix A. References 9 to 12 contain detailed discussions 
of the method of solution. All the measured quantities used in the anal- 
ysis. and the resulting aerodynamic derivatives are referred to the system 
of body axes illustrated in figure 5- 

Time histories of two typical oscillations for each of the three 
models are shown in figures 6 , 7> and 8. For convenience in illustrating 
the important features of the motion, the mean values of rolling velocity 
have been subtracted from the data and only the disturbances superimposed 
on these mean values are shown in figures 6 to 8. The essentially steady- 
state or mean values of rolling velocity were caused by the slight amount 
of twist in the model wings. The most obvious differences between the 
motions for the three models are the greatly reduced oscillatory rolling 
velocities for models 2 and 3 compared with those of model 1. At speeds 
below about M = 0 . 9 , model 3 exhibited no discernible periodic rolling 
motion. Also evident from figures 6, 7 , and 8 is the lower frequency of 
the motion for model 1 compared with that of models 2 and 3 * 

In order to illustrate the effects of these motion characteristics 
on the analysis of the data herein, vector diagrams for M 1.2 are 
shown in figure 9 for the high-wing model (model l) and the low-wing 
model (model 2) . The larger size of the side-force and yawing-moment 
diagrams for model 2 than for model 1 is a direct result of the higher 
frequency for model 2. The larger size of the rolling-moment diagram 
for model 1 is caused by its higher ratio of rolling to sideslipping 
motion. The side-force diagram shows that ft and -r are very nearly 
equal in magnitude and differ in phase angle by only a few degrees; 
these conditions justify the use of ft = -r in the rolling-moment and 
yawing -moment e quat i ons . 
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For solution of the rolling-moment equation, it is necessary to know 
one of the aerodynamic derivatives before the other two can be found. For 
the high-wing model, C i r - Cjjj was estimated and then and C 

were determined from the vector diagrams. It can be seen from figure 9(b) 
t ha t C i r - C has a relatively minor effect on the motion of the high- 

wing model. When the same procedure was used for the low-wing model, very 
great scatter in the values of Cjp occurred and some improbably high 
and some improbably low values resulted. It is evident in figure 9(1) 
t ha t C^ r - C^p has a much greater effect on the motion of the low -wing 

model. Comparisons of the two configurations (both at an angle of attack 
of nearly 0 °) indicated that the contribution of the wing-body combination 
to Cjp should be the same for the two models. Estimations of the con- 
tribution of the vertical tails to C^ (ref. 14) showed some differences 

for the two models because of the difference in sidewash produced by the 
rolling wings, but the entire contribution of the vertical tail was less 
than 5 percent of that due to the wing. Therefore, the values of Cjp 

for the two configurations should be the same within a few percent. On 
the basis of this reasoning the values of Cjp found for model 1 were 

assumed to apply for model 2 and the vector diagrams were completed to 
determine Ci r - Cjp and C^p. 

In the yawing-moment diagrams (fig. 9(c)), it can be seen that the 
product of inertia and Cn terms have an appreciable effect on the 

motion for model 1 , and a two-degree-of -freedom analysis of the yawing 
motion would be considerably in error. For model 2, however, neglecting 
the effects of the rolling motion would give a very good approximation 
to the yawing-moment derivatives. For solution of the yawing-moment 
equation, estimated values of = 0.02 for model 1 and C n p = 0.05 

for model 2 were used. These values were estimated by use of refer- 
ences 14 to 16 and the difference in value for the two models is due to 
the difference in sidewash at the vertical tail caused by the rolling 
wing. Reference 17, containing values of for a swept-wing con- 
figuration at subsonic speeds, showed values of no larger than 0.02 

from M = 0.4 to M = 0.95* 

For model 3 "the analysis was similar to that for model 2 except that 
below a Mach number of about 0.95 where the rolling motion became essen- 
tially zero, all the rolling terms disappeared from the equations and it 
was not necessary to estimate any of the aerodynamic derivatives. 

The measured values of Cy were plotted against 3 and the slopes 
of these curves yielded the derivative Cy^. Since the phase angle 
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between Cy and 0 was l80° for all models within the accuracy to whi ch 
this quantity could be measured, the amplitude ratio F— I is equal 


to 



The flexibility of the wings and vertical-tail surfaces of the models 
had an appreciable effect upon the results. The effects of wing flexi- 
bility on Cj and the increment in caused by geometric dihe dral 


and the effects of vertical tail flexibility on 
estimated by the methods outlined in appendix B. 



C np , and C nr were 


Accuracy 

An investigation of the accuracy of the lateral stability deriv- 
atives was made as in reference 9* Two Mach numbers. ( 1.3 and 0.7) were 
selected and the investigation was made for models 1 and 2. The results 
for model 2 should also apply approximately to model 3. 

From previous experience with similar instrumentation and by inspec- 
tion of the data described herein, the accuracies of the basic measured 
quantities were estimated and are shown in table IV. The increments in 
the various derivatives caused by errors of these magnitudes were then 
determined mostly graphically from the vector diagrams. The resulting 
increments are shown in table V. The probable errors in table V are 
given in terms of both absolute magnitude and as percentages of the deriv- 
atives. The percentage errors have little meaning in some cases, such as 
for Cjp of model 2 for example, where the absolute error is small but 

the percentage error is large because C^p is near zero. Also shown in 

the lower parts of table V are the increments caused by certain amounts 
of uncertainty in the aerodynamic derivatives that had to be estimated 
to permit determination of the other derivatives. 


From table V it is evident that the static sideslip derivatives are, 
in general, more accurately determined than are the rotary derivatives. 
When the effect of dynamic pressure inaccuracies which have an appreciable 
effect on all derivatives at low speeds is excluded, the values of c np 

are affected greatly by inaccuracies in Pj the values of Cn r - Cnp, by 


inaccuracies in T \J2} "the values of Cjp, by inaccuracies in and 

When the rolling 


the values of C^ r - "by inaccuracies in 




velocity is large, additional large effects on u n p 


caused hy inaccuracies in € and a, and Ci 


Cr 
and 


are 


°n r ' S 

are also affec- 


and 

■p Cz p 

ted greatly by inaccuracies in fip. An understanding of the reasons for 

all these effects can be obtained by a study of the vector diagrams in 
figure 9- 
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The effect of uncertainties in the estimated value of Cnp is fairly 
large for the yaw damping C nr - C n p for model 1 and fairly small for 

model 2, even though the estimated value is larger for model 2. This 
effect, of course, is again due to the great difference in the amplitude 
of the rolling motion between the two models. Errors in the values 
of C^p and C^ r - used in the analysis have appreciable effects 

on each other. However, because of the method of obtaining these deriv- 
atives for the present series of tests, as explained in the Analysis sec- 
tion, these quantities should be subject to less uncertainty than are the 
values of Cnp. 

RESULTS AND DISCUSSION 
Time Histories 

Portions of the model time histories were presented previously in 
figures 6 to 8 and seme characteristics of the lateral motions were 
pointed out. Small variations, mostly of a random character, occurred in 
the angle-of -attack and lift-coefficient time histories } these variations 
raise the question of inertia coupling between longitudinal and lateral 
motions. Although the effects of the lateral motion on the longitudinal 
motion were important in producing or modifying the longitudinal motion 
shown in figures 6 to 8, the factor of significance for the data pre- 
sented herein is the secondary effect of this longitudinal motion on the 
lateral motion. Calculations indicated that these latter effects were 
within the accuracy of the measurements for the data presented herein. 


Trim Characteristics 

The trim values or mean values of the measured quantities with oscil- 
latory components removed are shown in figure 10. These mean values are 
a measure of the asymmetry of the models or the measuring instruments. In 
particular, the positive value of Cy T for model 1 is believed to repre- 
sent a constant-acceleration error (or zero shift) in the transverse 
accelerometer at the center of gravity. This supposition is supported by 
the nonzero mean values which occurred in the values of i|r for this model, 
as calculated from the difference in readings between the accelerometers 
at the center of gravity and in the nose. Both of these discrepancies 
are of a magnitude and direction to be explained by the acceleration 
error mentioned. 

All models remained near zero lift. The mean values of rolling 
velocity are given in terms of the helix angle pb/2V, which was less 
than 0.01 for all models. The estimated a for model 1 shows the values 
used in the solution of the equations and was obtained by dividing the 
measured by an estimated Cjf a , since a was not measured on this 

model. 
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Oscillation Characteristics 

The measured periods of oscillation are shown in figure 11. The most 
obvious feature is the larger periods for the high-wing podel. The damping 
characteristics are shewn in figure 12. Little difference in magnitude 
existed between the damping times for the three models. 

The amplitude ratios of rolling velocity to sideslip angle are shown 
in figure 13 . The ratio for the high-wing model is 5 to 10 times as high 
as for the other two models. For both high-wing models, the rolling- 
velocity amplitude ratio increases considerably as the Mach number 
increases, whereas for the low-wing model it is more nearly constant. 

The 10° cathedral in the high-wing caused roughly the same decrease in 
rolling velocity as moving the wing to the low position. 

The phase angle by which the rolling velocity leads the sideslip 
angle is shown in figure 14. Below M = 0.90, no values could be deter- 
mined for model 3 because of the lack of any measurable oscillatory rolling 
velocity. Model 2 exhibits considerably larger variations of phase angle 
with Mach number than does model 1. At least part of this variation can 
be attributed to a lesser degree of experimental accuracy in determining 
the phase angles for model 2 because of the smaller amplitudes of rolling 
motion. This is indicated in table IV by the smaller accuracy listed for 
model 2. 

Plots of Cy against 3 are shown in figure 15 . Within the random 

scatter of the data the curves appear to be linear for models 1 and 2 . 

For model 3 some possible nonlinearity is evident in the ifa-ch number 
region from 0. 7^ to O. 83 , although linear slopes are faired through the 
data points in figure 15(c) for purposes of data analysis. The more acute 
angle at the wing-fuselage juncture for this model may be a source of flow 
separation leading to possible nonlinearities. No fairings of the wing- 
fuselage juncture were used on any of the models. As indicated by the 
vector diagrams in figure 9 , the effects of change in slope of the curve 
of Cy plotted against 3 will be very small. The main significance 
of any nonlinearities in the curve of Cy against 3 is that this result 
indicates the possibility of nonlinearities in other derivatives. Because 
the analysis used herein is a linear analysis, the derivatives obtained 
represent averages over the sideslip range of the tests. 

From the previous plots, slopes of the side-force coefficient as a 
function of sideslip angle have been determined and the derivative 
as a function of tfeich number is presented in figure 16. The high- 
wing model, as would be expected, has a lower Cy^ than does either of 

the other models. In the transonic region the low-wing model has a value 
of Cyp 25 to 50 percent greater than the high-wing model. The values 
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of Cyp for the high-wing model with cathedral fall between those for 

the other two models, these values being generally closer to those for 
the low-wing model. In the region of possible nonlinearity, the 
slopes Cyp for model 3 exhibit some rather large and irregular vari- 


ations with Mach number. As indicated in figure 1 6 , vertical- tail 
flexibility had a relatively small effect on Cy^. 


The preceding plots present the basic information which, with the 
mass and inertia characteristics given in tables I and II, constitute 
the data necessary to construct the vector diagrams. Quantities from 
the plots of period, amplitude ratio, and phase angle were not faired 
before using them to construct the vector diagrams since each point in 
effect constitutes a faired value to the extent that periods, amplitude 
ratios, and phase angles are for a number of half cycles of the oscil- 
lations. A faired value was used for T 1/2 in the data reduction because 
of the difficulty in obtaining an accurate value of this quantity. 


Figure 17 shows a comparison of the amplitude ratio of yawing accel- 
eration to sideslip angle for two models as measured from the transverse 
accelerometers by the relation 

r = C^nose a ^cg} g 

and as determined from the vector-diagram solutions. The agreement shown 
is considered to be very good and furnishes a check on the accuracy of 
the solutions. The comparison in figure 17 also is a further indication 
of the small magnitude of coupling effects because the difference in 
reading between two transverse accelerometers separated along the X-axis 
will actually be proportional to the sum (f* + pq) . 


Lateral Stability Derivatives 


The sideslip derivatives are presented for 0 in degrees to facil- 
itate comparisons with wind-tunnel test results. In the equation of motion 
these derivatives must be used with 0 in radians. The static stability 
derivative Cnp is presented as a function of ftkch number in figure 18 . 

In agreement with the high value of sidewash expected for a high-wing 
configuration (ref. l 8 ), the static stability derivative, of the high-wing 
model is consistently lower than that for either of the other two models. 

An approximately constant difference exists between the high- and low- 
wing models. The static stability derivative of model 3 is slightly lower 
than that for the low-wing configuration at subsonic Mach numbers but is 
approximately the same at supersonic Mach numbers. Hence, it appears that 
the lower values for Cy^ for model 3 than for model 2 (fig. l6) at super- 
sonic speeds may not be a result of sidewash effects on the vertical tail 
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but may "be a wing-fuselage effect. The effect of vertical-tail flexi- 
bility is greater for tbe low wing than for the high-wing model because 
at a given sideslip angle the vertical-tail load is greater for the low- 
wing model as a result of the smaller sidewash. 


The variation of the effective dihedral derivative C 7 „ with Mach 

tp 

number is presented in figure 19. The high-wing model had a large neg- 
ative value of which rose to a peak at M = 1.0 and decreased in 

magnitude at supersonic speeds. The values of Cjp for the low-wing 

model were approximately one -fourth as large as those for the high-wing 
model. The variation of the derivative with Mach number for the model 
with cathedral was such that the values of C^p were lower than that for 

the low-wing model at subsonic speeds but increased greatly at Mach 
numbers above 1.0 and became greater than that for the low-wing model at 
the highest speeds attained in the tests. 


An exact estimate of the effects of flexibility on 
involved since both wing and vertical tail contribute to 


would be very 
Cip, and the 


wing load distribution and resulting sidewash changes as affected by wing 
height are unknown. The increment in C contributed by the vertical 


tail will be reduced by vertical-tail flexibility but this reduction will 
be small (about 9 percent at the highest Mach number ) and will be nearly 
the same for all three models. The increment in caused by moving 

the wing to the top or bottom of the fuselage is essentially a wing-root 
effect and the influence coefficients in figure 3 indicate very little 
wing twist due to loads concentrated near the inboard end. Therefore, to 
a first approximation, the values of for the high- and low-wing 


models should be relatively unaffected by wing and tail flexibility. The 
increment in caused by the geometric dihedral is shown by refer- 

ence 19 and data to be presented subsequently herein to be approximately 
proportional to Cjp, although, of course, the load distributions are not 


identical. Therefore, as a gross approximation, the increments in C 


*3 


due to geometric dihedral were corrected by the same factor calculated 
for C 7 and the results are shown in figure 19 . 


As discussed in the Analysis section, the damping-in-roll deriv- 
ative C 7 was determined for model 1 only and it was assumed that 
4 P 

the C^p for models 2 and 3 was the same. Accordingly, the variation 
of Cj with Nfech number for the models is shown in figure 20. The 
effect of wing flexibility on Ct was very large, the estimated incre- 

lr 

ment due to elasticity being about 50 percent of the measured value 
at M = 1.3. The Cj data are almost entirely a wing-fuselage effect 
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since, as mentioned previously, the tail surfaces contribute less than 
5 percent of the damping in roll when the effect of the rolling wake on 
the tail is accounted for. 


Presented in figure 21 is the variation of the damping-in-yaw deriv- 
ative C nr - C n ^ with Mach number. As might be accounted for by 

the high rate of change of sidewash, the damping- in-yaw derivative for 
the high-wing model is generally higher than the damping for either of 
the other two models, although the differences approach the order of 
accuracy (table V) . Both models 2 and 3 show irregular variations in 
the supersonic region. No large loss in damping occurred near a Mach 
number of 1.0 as has generally been experienced for the pitch damping 
of sweptback configurations (ref. 20) . Flexibility of the vertical 
tail caused an estimated loss in yaw damping of about 9 percent at the 
highest Mach number. Since this value is less than the calculated 
probable accuracy of the derivatives, it is not shown in figure 21. 


The variation of the cross derivative C2 r - with Mach 

number is shown in figure 22. The solution for the model with cathedral 
at subsonic speed is relatively independent of the assumption for 

since the excitation of the roll approached zero below M = 1.0. 

Above M = 1.0, the Cj^ - for this model rises rapidly. For the 

low-wing model, the cross derivative has an erratic variation through 
the transonic region with a slightly negative value in the region 
around M = O.85. Because of the relatively large possible inaccuracies 
in this derivative (table V) and the unknown aerodynamic interference 
effects, no attempt was made to estimate the effects of flexibility. 


Comparison With Other Data 

Figures 23 to 25 contain summary plots comparing the data described 
herein to various wind-tunnel, rocket-propelled model, and theoretical 
results on similar configurations. The rocket-propelled model data of 
reference 21 are for wing and tail geometry identical to the models of 
the present tests except for a midwing position. The fuselage shape was 
somewhat different. The wind-tunnel data at Mach numbers of 2.01 
(ref. 22) and 1.4l (unpublished) were for the same wing as the present 
tests but a different empennage arrangement. (See fig. 23.) The wind- 
tunnel data at M = 0.17 (ref. 18) and M = 0.80 to 0-92 (ref. 23) 
were for configurations with somewhat different wing and tail geometry. 
All these wind-tunnel data have been corrected for the differences in 
tail geometry for comparison with the present tests. 

The curves shown in figures 23 to 25 for the present tests are those 
labeled "estimated rigid" in previous figures since the wind-tunnel data 
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are for essentially rigid models. In addition, the values of for 

the present tests have been adjusted to a center-of -gravity position 
at 0.20c to furnish a more valid comparison with data in reference 21 
(center of gravity at 0 . 19 c) and references 18 and 22 (center of gravity 
at 0.255). 


Results for the midwing rocket-propelled model of reference 21 agree 
very well with the sideslip derivatives of the present test (fig. 25 ) 
except for at M = 1.3 where the values from reference 21 are 

somewhat higher than would be expected for a midwing configuration from 
the present tests. The increment in and caused by changes in 

wing height, at an angle of attack of 0°, are caused by the vertical tail 
being in different regions of the sidewash field for the different wing 
positions. These increments for the present test are about the same as 
those obtained from the subsonic wind-tunnel data of reference 18 where 
a somewhat similar vertical tail was used and are somewhat larger than 
the increments from reference 25* The supersonic wind-tunnel data 
at M = 1.4l and M = 2.01 show much smaller increments resulting from 
changes in wing height. This condition is caused by the smaller changes in 
sidewash on the vertical tail for these tests as a result of the vertical 
tail being effectively nearer the center of the sidewash field. A vertical 
tail having its area symmetrically distributed above and below the longi- 
tudinal body axis would produce no change in Cy^ and C n g at an angle 

of attack of 0° as the wing is moved from the bottom to the top of the 
fuselage. 


The increments in 



caused by changes in wing height are prima- 


rily an effect of the wing and secondarily, an effect of the tail. These 
increments are reasonably constant throughout the Mach number range 
from 0 to 2.0, some increase occurring at transonic speeds. The comparison 
between the various sets of data for C^g in figure 23 is very good. No 

theoretical values are shown in figure 23 because no theoretical values 
of sidewash resulting from the load distributions due to wing-height 
changes are available. 


A summary of the effects of geometric dihedral on C^g is shown in 

figure 24. The approximate method of reference 19 when used with the 
experimental (flexible) values of CS from the present test showed 

Jr 

remarkable agreement with the experimental values. This fact formed the 
basis for the method of correcting the C^g value of the high wing 

with cathedral model for flexibility as described previously. The sub- 
sonic theory of reference 24 and the supersonic wind-tunnel data of 

c *p 

reference 25 show values of — — similar to those for the present tests. 
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The values of from the present test corrected to rigid-wing 

values (fig. 25 ) are somewhat higher at subsonic speeds than those 
attained from wind-tunnel data (ref. 17 ) or theoretical values (ref. 24), 
are considerably higher than the theory of reference 26 at M = 1.0, 
and are in very good agreement with the theoretical values in ref- 
erence 27 at supersonic speeds. The unpublished rocket-propelled-model 
data in figure 25 are for a midwing model with a wing having the same 
aspect ratio and sweepback as the wing of the present tests but less taper 
and thicker airfoil. These data show good agreement with the present 
tests at supersonic speeds but show a lower peak at M = 1 . 0 . The curve 
for the present test, when extrapolated, would indicate a slightly 
higher than the data in reference 21 which were for the same wing 

and tail arrangement. This is reasonable since the data of reference 21 
were not corrected for flexibility. 


Little information exists on the derivative C^ r - for com- 

parison with the present tests. Calculated values from reference 15 
at low subsonic speeds were about 0.04 which is at least the correct 
order of magnitude. (See fig. 22.) Reference 6 contains measurements 
of C^ r - at supersonic speeds for a swept-wing configuration, and 

the values obtained were much larger than those estimated, the reason for 
this not being obvious from the tests. 


Flying Qualities 


The motion characteristics presented previously for the three models 
may be interpreted in terms of full-scale airplane flying qualities. The 
nondimens ional characteristics of the motion, that is, amplitude ratios 


such as I 


4 


phase angles, and number of cycles required to damp to a 


given fraction of the original amplitude may be considered as applying 
directly to a full-scale airplane having the same relative density and 
the same nondimens ional radii of gyration. 


The relative density factor for the model tests varied from 

about 115 at a Ivfe.ch number of 0.7 to about 87 at a Mach number of 1 . 3 . 
The full-scale airplane altitudes and wing loadings represented by these 
values, the models flown being assumed to be l/lO-scale models, are 
shown in figure 26. The values simulated by the models are reasonable 
for current aircraft. 


Requirements for damping of the lateral oscillations for military 
aircraft are given in reference 28. These requirements are illustrated 
in figure 27 . Points for all three models discussed herein are shown 
on the chart for two Mach numbers. It is apparent that the damping is 
marginal or unsatisfactory for all three configurations for normal-flight 
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conditions and would require augmentation to satisfy tactical-mission 
requirements (gunnery or bombing runs). 


A comparison was made of the nondimens ional radii of gyration of the 
three models with average values for three current swept-wing fighter 
aircraft. The comparison showed that the radii of gyration in roll were 
the same but that the radii of gyration in yaw of the models were about 
15 percent larger than those for the full-scale airplanes. Also, as 
stated previously, the centers of gravity of the models were placed near 
the leading edge of the mean aerodynamic chords to minimize pitch-yaw 
coupling. A normal center-of -gravity position for this type of aircraft 
would probably be near 30 percent of the mean aerodynamic chord. The 
effects on the motion of reducing the radius of gyration in yaw by 15 per- 
cent and the effects of this reduction combined with a 30 -percent rearward 
movement of the center of gravity were investigated for the high-wing 
model at a Mach number of 1 . 3 . This point was chosen since it represents 
the lowest static stability measured during the present tests and thus 
would be expected to produce the largest effects. 'Hie stability deriv- 
ative C n p and the damping derivative C nr - C^. were transferred 

to the altered center-of -gravity position. Since the transfer equations 
for C nr and Cnjj are different, this procedure cannot be done rig- 


orously since only the quantity C nr - is known.. Wind-tunnel tests 

and est ima ted results, however, indicate that should be small com- 

pared with C nr at low lift and that C nr is nearly all caused by 
the vertical tail. Therefore, the damping derivative C nr - Cn* was 
reduced by the square of the ratio of vertieal-tail distances from the 
two center-of -gravity positions. Solutions for these altered conditions 
were made with the aerodynamic derivatives as known quantities to obtain 
the motion characteristics. The results are shown in figure 2J. Reducing 
the radius of gyration caused the motion to move further into the satis- 
factory region. When this condition was combined with a rearward center- 
of -gravity movement, the motion became definitely unsatisfactory. This 
latter result is contrary to what might be expected if a simple yawing 
oscillation is considered. The calculations showed that it was caused 
by the product-of -inertia effect assuming relatively greater importance 
when the aerodynamic derivatives and C nr - C n £ were reduced by 

the rearward center-of -gravity movement. 


An additional desirable quality for the lateral oscillatory motion 


is that the ratio 


should not exceed about (ref. 29 ) to preclude 


excessive rolling motions resulting from small corrections in heading 
or from flight through rough air. Figure 28 shows the values of 
the three models. For the high-wing model, the ratio |— | remains 


IP 


for 
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above 5 for all ftkch numbers. For the other two models, 
below 2 over the Nfe,ch number range Investigated. 


t 

3 


remains 


CONCLUSIONS 


Lateral stability investigations were made at low angles of attack 
by flight tests between Mach numbers of about 0.7 and 1.3 of three rocket- 
propelled models of a swept-wing airplane configuration, one having a 
high wing, one a low wing, and one a high wing with 10 ° cathedral. 

Analysis of the motions following yawing disturbances and comparisons 
with other investigations indicated the following conclusions: 

1. Increments in the static lateral stability derivatives caused 
by moving the wing from a high to a low position were roughly constant 
in the transonic region and were of the same magnitude as the increments 
estimated for these configurations based on low-speed wind-tunnel data 
for a similar configuration. 

2. Estimates based on supersonic wind-tunnel data showed smaller 
increments due to wing-height changes on the static lateral-force and 
yawing -moment derivatives than the present tests. This difference 

can be explained by the differences in vertical-tail geometry resulting 
in different magnitudes of sidewash effect. 

3 . The change in the static rolling -moment derivative due to wing 
height, which is primarily a wing-fuselage effect, was of a similar 
magnitude for the present tests and the supersonic wind-tunnel tests 
which used the same wing and a similar fuselage. 

b. The effect of geometric dihedral on the static rolling -moment 
derivative showed good agreement with subsonic theoretical results and 
with a theoretical relationship between this derivative and the exper- 
imental damping in roll. 

5. The damping in roll, when corrected to rigid-wing values and 
compared with theoretical and other experimental data (on midwing config- 
urations) showed very good agreement at supersonic speeds but was some- 
what high at subsonic and transonic speeds. 

6 . When interpreted in terms of full-scale airplane flying qualities, 
the damping of the lateral oscillations of these configurations was not 
sufficient and augmented damping would be required. In addition, the 
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ratio of rolling motion to sideslipping motion for the high-wing config- 
uration was undesirably large. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , May 4, 1956. 
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APPENDIX A 
EQUATIONS OF MOTION 


All the measured quantities used in the lateral-stability analysis 
herein and the resulting aerodynamic derivatives are referred to the 
system of body axes illustrated in figure 5* 'Hie three equations of 
motion used to represent the free-oscillation response of the models in 
the lateral mode following a disturbance were: 


mV 
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X x 

qSb 


p _ r _ C 7 0 - Ct 2k _ fc _ C-,.\ 
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! lr “ C *p) 
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2V 

(A2) 

' n r " C np) 

— = 0 
2V 

(A3) 


In the side-force equation (Al) the gravity terms (0 cos 0 4 - sin 0) 

qb 

have been omitted. As shown in the Accuracy section of this report, this 
omission has a negligible effect on the results presented herein and 
eliminates the necessity of measuring or computing the roll, pitch, and 
yaw angles. Also in equation (Al) the total aerodynamic side force is 

given by Cy an( i includes the contributions of Cy_, Cy* , Cy , 

1 A p p p 

and Cy because this is the quantity obtained from the transverse 

r W 

accelerometer by the relation Cy = a t — . 

In equations (A2) and (A3), the assumption has been made that p =*-r 
in order that the yawing- and sideslipping-velocity derivatives may be 
combined to reduce the number of unknown aerodynamic terms. As shown by 
the vector diagrams in figure 9* this is a very good approximation for the 
tests discussed herein. 


Detailed discussions of the application of the time-vector method to 
the solution of equations (Al) to (A3) niay be found in references 9 to 12. 

The vector solution of the equations of motion is illustrated in 
figure 9. All quantities are referenced to p having an amplitude scale 
shown on the diagrams and zero phase angle. Only two of the three deriv- 
atives in each of the rolling-moment and yawing-moment equations can be 
determined from one set of test conditions by the method of solution 
described herein. The methods of approximation used in the solution are 
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discussed in the Analysis section of this report and some comments on 
extension of the test technique to include other modes of motion and 
other test conditions are included in appendix C. 

For model 3 the oscillatory rolling velocity was zero at subsonic 
speeds and for this case equations (Al) to (A 3 ) simplify to 

3T (5 + r) - C Y = 0 (A 4 ) 


T 

^r-Cn p 3-(Cn r -C n ^)-=0 


(A 5 ) 

(a6) 


Since Cn 


and 


Cn are eliminated from the equations, it was not 

.b' 


necessary to estimate any of the derivatives for solution of equa- 
tions (Ak) to (A 6) . 
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APPENDIX B • 

EFFECTS OF FLEXIBILITY OF VERTICAL TAIL AND WING 
Vertical Tail 


The effect of vertical- tail flexibility on the vertical- tail lift- 
curve slope was calculated by the method of reference 30. For this 
calculation the influence coefficients in figure 3 a nd the span- load 
distributions of reference 31 were used. The loss in vertical-tail lift- 
curve slope due to flexibility was 9 percent at a Mach number of 1.3 
and 1 percent at a Mach number of 0.7* Estimated rigid values for the 
derivatives Cyp and C n were then obtained from 



estimated from wind-tunnel tests. 


Wing 


The effect of wing flexibility on 


C 7 was calculated by the method 
ip 


of reference 30 with some modifications to obtain rolling-moment effects 
rather than lift effects. The influence coefficients used are shewn in 
figure 3 and the span- load distribution was obtained from reference 24. 


The equations herein represent the right wing of the airplane with 
a positive load L producing negative rolling moment. The wing and span- 
load distribution are divided into a number of sections, as in reference 30. 
The for the elastic wing is given by 
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The rolling moment produced by section 1 of the elastic wing can be 
approximated by 


*1 

where 


r( P ) 

span-load 

Axi 

change in 

’ll 

value of 
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From equation (B2) 



parameter for rolling wing 

section angle of attack caused by wing elasticity 
T] at centroid of load distribution over section 1 
r) at root of exposed wing 
T) at inboard end of section 1 
T) at outboard end of section 1 
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From the influence coefficients the change in angle of attack at sec- 
tion 1 may be obtained as 


/Vtn Lt Lo 

1 — = 1 — 6 11 ■* — ®12 + 

pb/2V pb/2V pb/2V 


(b6) 


where 9^ is the twist at station i due to a unit load at station j. 
Equating (B5) and (b 6) and rearranging yields 
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Similar equations are written for each section and the resulting set of 


equations is solved simultaneously for the values of 


Ll 


L 2 


pb/2V' pb/2V 


, and 


so forth, for various values of r q covering the test range. Then 


L 1 


and so 


from equations (Bl) and (b 7) and the values of , , 

pb/2V pb/2V 

forth, the ratio of elastic damping in roll to rigid damping in roll is 
given by 


Cl. 
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V,r 


I Iq L g 

\pb/2V 111 + pb/2V n? ~ 


+ . . 


Ci qS 
ip,r H 


(b8) 


In addition to its use in correcting Ci for wing flexibility, 

ir 

equation (b 8) was used to give a first approximation to the effects of 
wing flexibility on the increment of C ^ caused by geometric dihedral* 

as explained in the text. 
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APPENDIX C 

EVALUATION OF TEST TECHNIQUE 


Analysis of the data for the three configurations discussed herein 
provided some comparative hasis for evaluating the test technique. The 
most obvious conclusion is that no general statements applicable to all 
cases can be made and each case must be examined individually to determine 
the aerodynamic derivatives that can be determined and their resultant 
accuracy. 


The primary weakness in the test technique is that not all the 
rolling-moment and yawing-moment derivatives usually included in the 
equations of motion can be determined. This is a result of the fact 
that the Dutch roll lateral oscillation when analyzed by the time- 
vector method does not by itself furnish sufficient information for 
determining all the derivatives. Some approximations must be made to 
permit the determination of the most important derivatives. One approx- 
imation generally made is that p = -r which reduces the number of 
unknown derivatives by 2. 


In many cases the cross derivatives C nr) and C 7 - C 




have 


-np ^l r 

relatively small effects on the motion and their magnitude may safely be 
estimated to permit. determination of the other derivatives. For some 
cases where the effects of Cn or - Cj* may be appreciable but not 


necessarily larger than the effects of other derivatives, it may be advan- 
tageous to use estimated values of some of the other derivatives such 
as Cjp, Cjp, or Cn. r " Cnjjj because of the relatively greater amount of 

information available on these derivatives. 


It is also fairly obvious that the rolling derivatives C 7 and C 7 

L V L [ 

are best obtained from motions involving a large ratio of rolling to 
yawing motion and the yawing derivatives and Cn r - C n £ are best 

obtained from motions having a large ratio of yawing motions to rolling 
motion, because these conditions tend to minimize the effects of the 
cross derivatives. Conversely, the effects of the cross derivatives can 
be increased by reversing these conditions. Thus, intelligent planning 
of the dynamic conditions of the tests can result in emphasizing the 
derivatives of most interest. 


Alterations in the ratio of rolling motion to yawing motion may be 
accomplished by changing the ratio of rolling moments of inertia to yawing 
moments of inertia. If two flights at two different inertia ratios but 
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under otherwise identical conditions can be made, then the difficulty 
of having too many derivatives to measure can be obviated, because two 
different sets of motions would be available. Both sets of motions must 
be produced by the same se;t of derivatives, and thus a larger number of 
derivatives may be determined. 

In order to assist in analyzing the data other modes of motion 
may be used, such as the heavily damped roll subsidence. The damping 
of this motion is almost completely dependent upon C^, the other 

derivatives having only a negligible effect, and may be used to 

find C 7 if this mode of motion can be extracted from the measured 
P 

time histories. This procedure has been used in other unpublished 
rocket-propelled model tests to measure . For the present tests 

the time to damp to half* amplitude of this mode of motion was about 
0.03 second at M = 1.3 and 0.06 second at M = 0.7 and examination of 
the time histories indicated that this motion could not be extracted from 
the results. Greater excitation of this mode of motion can be obtained by 
applying a rolling-moment disturbance rather than a yawing-moment dis- 
turbance as for the present tests. 

Although the test technique used is a dynamic technique, the analysis 
of these and other flight test data indicated that the so-called static 
sideslip derivatives can usually be determined more accurately than can 
the rotary derivatives. The rotary derivatives depend primarily on the 
phase relations between the various components of motion and a more 
accurate determination of the rotary derivatives would require greater 
accuracy in the measurement of the phase angles and of those quantities 
which effectively cause phase shifts in the motion, such as angle of 
attack, inclination of the principal axis, and the damping of the motion. 
The requirement for accuracy of measurement of the phase angles also 
requires that the frequency-response characteristics of each instrument 
be known, since appreciable errors in phase angles due to imperfect 
instrument response may be present in the data and must be corrected for. 
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TABLE I 

GEOMETRIC CHARACTERISTICS OF THE MODELS 

Wing: 

Area (extended to model center line), sq ft 4.4l 

Span, ft A. 20 

Aspect ratio 4.0 

Mean aerodynamic chord, ft 1.21 

Sweepback ( quarter chord) , deg 45 

Taper ratio 0.20 

Airfoil section NACA 65A004 

Horizontal tail: 

Area (extended to model center line), sq ft 0.88 

Aspect ratio 4.0 

Span, ft 1.88 

Mean aerodynamic chord, ft 0.54 

Sweepback ( quarter chord) , deg 45 

Taper ratio 0.20 

Airfoil section NACA 65A004 

Vertical tail: 

Area (extended to model center line), sq ft O .89 

Aspect ratio 1.59 

Height (above center line of model), ft 1.19 

Sweepback ( quarter chord) , deg 45 

Taper ratio 0.17 

Airfoil section NACA 65A004 
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TABLE II 

MASS AND INERTIA CHARACTERISTICS 



Model 1 

Model 2 


Weight, lb 

103.0 

111.0 

112.7 

Center -of -gravity position, percent c • • 

0 

3.8 

-0.6 

Moments of inertia: 




Pitch, slug-ft^ 

5.V7 

5.88 

5.90 

Yaw, slug-ft^ 

5.81 

6.14 

6.38 

Roll, slug -ft ^ 

0.6l 

0.66 

0.69 

Inclination of principal axis, deg ... 

2.15 

0.75 

1.60 

Product of inertia 

0.20 

0.07 

0.16 






TABLE III 


FUSELAGE ORDINATES 


Station, 

in. 

Diameter, 

in. 

0 

0 

3 

1.60 

6 

3.00 

9 

4.24 

12 

5.28 

15 

6.14 

18 

6.84 

21 

7.34 

24 

7.66 

27.8 

7.80 

30 

7.78 

33 

7.74 

36 

7.64 

39 

7.48 

42 

7.30 

45 

7.06 

48 

6.78 

51 

6.44 

54 

6.08 

57 

5.66 

60 

5.18 

63 

4.68 

66 

4.12 

69.5 

3.42 
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TABLE IV 


ESTIMATED ACCURACY OF VARIOUS MEASURED QUANTITIES 


^All increments may "be positive or negative] 


Quantity 



Accuracy at - 


M = 1.3 


M = 0.7 


M, percent 
q, percent 

Ti/2^ sec 


deg 


W, percent 

O 

k z , percent 
,, 2 


1 and 2 
1 and 2 
1 and 2 
1 and 2 


percent 


1 and 2 
1 and 2 

1 and 2 
1 and 2 
1 and 2 
1 and 2 

1 and 2 


1.0 

2.0 

.03 

.005 

10 

4 


3.0 

6.0 
.10 
.010 

3 

2 
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TABLE V.- CONCLUDED 

CALCULATED ACCURACY OE AERODYNAMIC DERIVATIVES 
(b) Model 2 


Increment due to error in - 

ACy 

3 

AC n 

P 

ac 2 

p 

A(C lr 

' Ci p) 

A(Cn r 

Cn &) 

M = 1.5 

M = 0.7 

M = 1.3 

M = 0.7 

M = 1.3 

M = 0.7 

M = 1.3 

M = 0.7 

M = 1.3 

M = 0.7 

W 

0.004 

0.004 

0.001 

0.001 

0 

0 

0 

0 

0.002 

0.003 

k 2 
K z 



.003 

.004 





.006 

.008 

v 2 





0 

0 

, oo^ 

.003 






0 

0 

.001 

.002 

.00^ 

.00U 

0 

.016 




0 

0 

0 

0 

0 

0 

.020 

.016 

at 

B 

.018 

.014 

0 

0 

0 

0 

0 

0 

0 

0 


.018 

.043 

.004 

.014 

0 

.001 

.002 

.006 

.008 

.032 

Tl /2 



0 

0 

0 

0 

.001 

.002 

.063 

.093 

P 



.015 

.015 

.001 

.001 

.004 

.005 

.020 

.039 

P 

— 


0 

0 

.002 

.003 

.020 

.019 

0 

0 




0 

0 

.001 

0 

. 03 ^ 

.036 

0 

.010 

Omitting veight effect 



0 

0 

0 

0 

0 

0 

0 

0 

Provable error, \jz( A ) 2 

.026 

.045 

.016 

.021 

.003 

.004 

.o4o 

.042 

.070 

.108 

Value of derivative 

-.877 

-.716 

.212 

.241 

-.017 

-.020 

.120 

.100 

-.420 

-.540 

Error, percent of 











derivative 

3 

6 

8 

9 

18 

20 

33 

k2 

17 

20 

c Zp (0.05 error) 






.001 

.020 

.020 






.001 

• 001 





.013 

.020 

Lj^p ^ luo percent error / 
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Model 3 

(to) Front views of models. 
Figure 1.- Concluded. 
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Model I 




( a ) Sideviews . 

Figure 2.- Photographs of the models. 
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(b) Plan view of model 2. ^“93515 


Figure 2.- Continued 
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(c) Model 1 on the launcher. 
Figure 2.- Concluded. 
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(a) Wing. (b) Vertical tail. 

Figure 3*- Measured influence coefficients. 
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(a) Dynamic pressure. 



M 

(b) Reynolds number based on mean aerodynamic chord. 
Figure Ij-.- Test conditions. 
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(a) Supersonic oscillation. (b) Subsonic oscillation. 


Figure 6.- Time history of model 1 
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(a) Supersonic oscillation. 

Figure 7.- Time hi' 



of model 2 
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Figure 8 .- Time history of model 3 
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(a) Side -force diagram. 


Figure 9*- Vector diagrams for high-wing and low -wing models. M ~ 1.2. 
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(a) Model 1. High-wing model. 
Figure 10.- Model trim characteristics 
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(t>) Model 2. Low-wing model. 


Figure 10.- Continued. 
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Figure 11.- Period of the oscillation 
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(a) Model 1. High wing. 



(b) Model 2. Low -wing model. 



(c ) Model 3 . High-wing model, 10° cathedral. 
Figure 12.- Time to damp to one -half amplitude. 
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M=0.94 1.01 l.il 1.21 1.31 

(a) Model 1. High-wing model. 


Figure 15*- Lateral-force coefficients 
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(b) Model 2. 
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Low -wing model. 
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Figure 15.- Continued. 
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(c) Model 5 . High-wing models 10° cathedral. 


Figure 15.- Concluded. 







Figure 19- - Effective dihedral derivative 
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Figure 20.- Damping -in-roll derivative. 
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Figure 21.- Damping-in-yaw derivative. 
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Ref. 25 


Figure 2k. - Effect of geometric dihedral on the effective dihedral 

derivative C ^ . 
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Figure 25.- 


Comparison of rocket -propelled model, wind-tunnel, and 
theoretical damping in roll. 
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Figure 28.- Roll-to-sideslip ratios. 
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